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We theoretically study the spin-dependent transport in a ferromagnet/super- conduc- 
tor /ferromagnet double tunnel junction. The tunneling current in the antiferromagnetic alignment 
of the magnetizations gives rise to a spin imbalance in the superconductor. The resulting nonequi- 
librium spin density strongly suppresses the superconductivity with increase of bias voltage and 
destroys it at a critical voltage V c - The results provide a new method not only for measuring the 
spin polarization of ferromagnets but also for controlling superconductivity and tunnel magnetore- 
sistance (TMR) by applying the bias voltage. 
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Since the early experiments demonstrated the spin- 
polarized tunneling of electrons from ferromagnetic met- 
als (FM) into superconductors (SC) in FM/SC junctions 
jjj , the concept of the spin-polarized transport has been 
of vital importance in magnetic junctions and multilay- 
ers. Firstly the tunneling currents strongly depends on 
the relative orientation of magnetizations in FM/FM tun- 
nel junctions; the tunnel resistance decreases when the 
magnetizations are aligned in a magnetic field, causing 
tunnel magnetoresistance (TMR) [Eh5| . Secondly the 
spin-polarized current driven from a FM into a normal 
metal (N) or superconductor (SC) gives rise to a nonequi- 
librium spin density in N or SC §-|]. In a FM/N/FM 
double junction the TMR effect is brought about by ac- 
cumulation of spin-polarized electrons in N 1. In a 
FM/SC/FM double junction ||, on the other hand, we 
expect the strong competition between superconductiv- 
ity and magnetism induced by the spin polarization in 
SC. Of particular interest in the FM/SC/FM structure 
is not only to find novel magnetoresistive effects due to 
the competition but also application to magnetoelectron- 
ics. 

In this Letter we show that a FM/SC/FM double tun- 
nel junction is a new magnetoresistive device to control 
superconductivity by applying the bias voltage (or cur- 
rent). In the antiferromagnetic (A) alignment of mag- 
netizations (see Fig. 1), a nonequilibrium spin density is 
induced in SC due to the imbalance in the tunneling cur- 
rents carried by the spin-up and spin-down electrons, so 
that the superconducting gap A is reduced with increas- 
ing bias voltage and vanishes at a critical voltage Ve- 
in the ferromagnetic (F) alignment, however, there is no 
spin-density in SC. Consequently, TMR has a strong volt- 
age dependence around V c ; TMR is enhanced compared 
with that in the normal state above V c , while it changes 
sign to show an inverse TMR effect for some voltage range 
below V c . It is shown that V c is inversely proportional 
to the spin polarization P of FM (V c oc 1/-P), which pro- 
vides a new method for determining P of FM. 

We consider a FM/SC/FM double tunnel junction as 



shown in Fig. 1. The left and right electrodes are made 
of the same FM and the central one is a thin film SC. 
The magnetization of the left FM is chosen to point up 
and that of the right FM is either up or down. The volt- 
ages —V\ and V% (= V — V\) are applied to the left and 
right electrodes, respectively. We assume that the en- 
ergy relaxation time te of quasiparticles in SC is shorter 
than the time Tt between two successive tunneling events, 
whereas the spin relaxation time r s is longer than r t . 
Consequently, the electrons tunneling into SC relax to 
the Fermi distribution before leaving SC, keeping their 
spin direction during the stay in SC. 
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FIG. 1. Double tunnel junction consisting of two ferro- 
magnets (FM) and a superconductor (SC) separated by in- 
sulating barriers in the case of identical barriers. Schematic 
densities of states of FMs (left and light) and SC (middle) 
in the antiparallel alignment of magnetizations in FMs are 
shown when SC is in the normal state (a) and in the super- 
conducting state (b). Sfi denotes the shift of the chemical 
potentials in SC. 
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We calculate the tunneling current using a phenomeno- 
logical tunneling Hamiltonian. If SC is in the super- 
conducting state, it is convenient to rewrite the electron 
operators in SC in terms of quasiparticle operators 
7ko- appropriate to the superconducting states, using the 
Bogoliubov transformation, a k ^ — Uk7kt + w k7-kj. ano - 
ak| = M-k7ki - v*_ k ll kV where |w k | 2 = \ (1 + Ck/^k), 
W\ 2 = |(l-W^k), and E k = [£ + A 2 ]V2 is the 
quasiparticle dispersion of SC, £k being the one-electron 
energy relative to the chemical potential and A being the 
gap parameter. Then, using the golden rule formula, we 
calculate the spin-dependent currents Ij a across the jth 
junction. The results are 



lit = (Gi T /e) [A/i - S - Q*/2] 
I n = (G n /e)[Af 1 +S-Q*/2] 
/ 2T = (G 2T /e)[7V 2 +5 + Q72] 
I 2l = (G 2l /e)[Af 2 -S + Q*/2] 



(la) 
(lb) 
(1c) 
(Id) 



where Gj a is the tunnel conductance of the jth junction 
for electrons with spin a if SC is in the normal state. The 
quantity Afj is given by the usual expression 1 1 1 



Nj = / V s (E k ) f (E k - eVj ) - f (E k + eV 3 ] 



dE k 



(2) 



where V s (E k ) = E k j ^[E\ - A 2 is the normalized BCS 
density of states and fo(E k ± eVj) is the Fermi distribu- 
tion function of thermal equilibrium in FM. The quantity 
S represents the spin density normalized to the density 
of states of SC in the normal state 2?n and is given by 



V s (E k )(fo-f kl )dE k , 



(3) 



where f kr7 — (7 ko .7 kcr ) is the distribution function of 
quasiparticles with energy E k and spin a in SC. The 
quantity Q* is the charge density normalized to eT>^ due 
to the imbalance in populations of electronlike and hole- 
like quasiparticles 



n_pc 



= 2 



(/kWk.) dE k 



(4) 



where f ka and f ka represent f kr7 in the electronlike (k > 
kp) and holelike (k < kp) branches of E k , respectively. 

In the limit of vanishing spin-flip scattering, the spin- 
up and spin-down currents are treated as independent 
channels. The conservation of the currents at junctions 
1 and 2, I\ a = /20-, yields the relations 

S = [(G l1 G 2l -G n G 21 )/(G 1 G L )] (A/i+A/- 2 )/2, (5) 
Q * = £ [( G i<r/Ga)Mi - [G 2a /G a )N 2 ] , (6) 



where G a = G\ a + G 2a . If the tunnel barriers are sym- 
metric [see Fig. 1], V x = V 2 = V/2 and M = M 2 = Af, 
so that Q* = for both alignments. This is because 
the charge transport is symmetric at the two junctions 
where the injected and extracted charges are balanced. 
In the following we assume the identical tunnel barriers 
and neglect the effect of the charge imbalance. 

The spin density S depends strongly on whether the 
magnetizations in FMs are parallel or antiparallel. The 
spin density Sa in the A-alignment (G\ a = G 2 - a ) satis- 
fies the relation 



S A = iW, 



(7) 



where P — \Gj-\ — GjA/ {G^ + Gj±) represents the spin 
polarization of FM [y. The relation (0) implies that a 
finite Sa is induced by applying the bias voltage. The 
Sa 7^ is a consequence of symmetry breaking of the spin 
transport in the A-alignment, and is realized when the 
distribution of the spin-up and spin-down quasiparticles 
is disequilibrium, i.e., f k ^ ^ f k ±, as seen from Eq. (|^). 
However, the spin density Sf in the F-alignment (Gi a = 
G 2a ) has no net spin-density in SC [Sp = 0) due to the 
symmetric spin transport at the junctions. 

The distribution function f ka is determined as fol- 
lows. When the thickness of SC is much smaller than the 
spin diffusion length |J, the distribution of spin- up and 
spin-down quasiparticles is spatially uniform in SC. For 
te < 7t < t s , the distribution is described by the Fermi 
function /o, but the chemical potentials of the spin-up 
and spin-down quasiparticles are shifted oppositely by 
Sfj, from the equilibrium one (see Fig. 1) to produce the 
nonequilibrium spin density. Thus we write f ka as ||ll[j 



/kt =fo(Ek-Sfx), f kl = f {E k + 8n). 



(8) 



In the normal state (A = 0), from Eqs. (|J), (0), and 
(H), we have 8^a — Sa — \PeV in the A-alignment [||, 
whereas 8\xf = Sf = in the F-alignment. 

We first discuss how the superconductivity is affected 
by the nonequilibrium spin imbalance in SC. The gap 
A in the nonequilibrium situation is determined by f ka 
through the BCS gap equation |To[ ] 



nVbcs 



1 - /kt - /kj 
E k 



(9) 



where f k(T is given by Eq. (§). We note that Eq. (|) with 
Eq. (||) is the same as that of SC in the paramagnetic 
limit if 5 /j, is taken to be the Zeeman energy hbH [ p^| - 
The chemical potential difference 25/i plays the role of 
pair breaking energy. Therefore the superconductivity is 
destroyed in the A-alignment when Sfj, exceeds a certain 
critical value by increasing the voltage. To show this, we 
solve self-consistently Eqs. ^ and (£)]) with respect to A 
and Sfi, and obtain A and S as functions of V. 

Figure 2(a) shows the gap parameter A^ in the A- 
alignment as a function of bias voltage V for P = 0.4, 
the spin polarization of Fe fl . The quantity A denotes 
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the value of for P = at T = 0. The gap parame- 
ter A,4 decreases with increasing V and vanishes at the 
critical voltage V c . At very low temperatures A a be- 
comes multi-valued in a certain range of eV just below 
2Ao; At T = it has t hree solutions, A^ = Ao and 
A A = A [l - 2P 2 ± 2P v / (eV/2A ) 2 + P 2 - l] 1 / 2 , in the 
range 0.92 < eV/2A n < 1. When V is increased (or 
decreased) at T ~ 0, an instability into a spatially inho- 
mogeneous state with different Aa takes place at a cer- 
tain voltage within the range. At V c where Aa — and 
<5/M — Eq. (^|) reduces to an implicit equation 

for PeV c and T, which gives a universal relation between 
PeVc/ 'Ao and T/T c , as shown in the inset of Fig. 2(a). 
In particular, we have eV c = A /P at T = 0. Therefore, 
we can determine the spin polarization of FM by mea- 
suring V c - Since the paramagnetic effect caused by spin 
accumulation becomes stronger with decreasing T, V c is 
not a monotonic function of T, but has a maximum at 
T/T c — 0.5. Figure 2(b) shows the voltage dependence of 
the spin density Sa in the A-alignment. The dotted line 
indicates the values of S A = \PeV in the normal state. 
As T is lowered below T c , Sa is suppressed below V c by 
the opening of the energy gap. At and near T = 0, Sa 
shows the S-shaped anomaly around eV c ~ 2Ao, which 
stems from the multiplicity of A^ shown in Fig. 2(a). 
The detailed behavior of the anomaly is shown in the 
inset of Fig. 2(b). In the F- alignment, Af has no V 
dependence and has the same value as Aa(V = 0). 

We now calculate the tunneling current as a function 
of bias voltage V. From Eqs. (2a)-(2d), the total currents 
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FIG. 2. (a) Gap parameter Aa in SC as a function of bias 
voltage V for different temperatures below T c in the antifer- 
romagnetic alignment. The inset shows the critical voltage 
Vc vs temperature T. (b) Spin density Sa in SC as a func- 
tion of V. The inset shows an enlarged view of Sa around 
eV/2A = 1 at T/T c = 0, 0.05, and 0.1. 



Jp and Ia in the F and A alignments are given by 

I F (V) 
Ia(V) 



(Gfn / e)Af(V, Ap), 
(G m /e)(l-P 2 )M(V,A A ), 



(10) 
(11) 



) and Gfn = Gj-\ + Gj±. 
) that I a oc Sa, so that 



where Af(V, A) is given in Eq. 
It follows from Eqs. (^) and flllQ 
Fig. 2(b) also represent the V dependence of 7a- 

Figure 3(a) shows the voltage dependence of the dif- 
ferential conductance Gf and Ga for the F and A align- 
ments at T/T c = 0.4. The Gp shows the ordinary de- 
pendence on V expected for the constant gap Ap- In 
contrast, because of the decrease in A^ with increas- 
ing voltage, Ga increases with voltage more rapidly than 
Gf, forming a higher peak than Gf, and then decreases 
steeply. At V c , Ga jumps to the conductance G A in the 
normal state. The tunnel magnetoresistance (TMR) is 
calculated by the formula: TMR = (G F /G A ) - 1- Using 
the values of Fig. 3(a), we obtain the V dependence of 
TMR shown in Fig. 3(b). At V = where A A = A F , 
TMR takes the same value as in the normal state. A deep 
negative dip appears at eV/2Ao ~ 1 where Aa steeply 
decreases, exhibiting an inverse TMR effect (Ga > Gf), 
and is followed by the discontinuous jump at V c above 
which TMR is highly enhanced compared to that in the 
normal state. 

The relation Ia oc Sa in the A-alignment directly in- 
dicates that the superconductivity of SC is strongly sup- 
pressed with increase of injection current I a- Using the 
relation and the result of A^ and Sa in Fig. 2, we obtain 
A,4 as a function of injection current I a- Figure 4 shows 
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FIG. 3. (a) Tunnel conductance as a function of bias volt- 
age. The dashed and solid curves indicate the conductances 
Gf and Ga for the ferromagnetic and antiferromagnetic align- 
ments, respectively, (b) Tunnel magnetoresistance (TMR) as 
a function of bias voltage. The dotted line indicates TMR 
= P 2 /(l — P 2 ) in the normal state. 
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FIG. 4. Square of Aa vs injection current I a at different 
temperatures below T c . A A is proportional to the superfluid 
density and thus represents the critical current of SC. 

the square of A^ as a function of I a- Since A 2 A is propor- 
tional to the superfluid density which is a measure of the 
critical current I c of SC, Fig. 4 represents also the critical 
current as a function of the spin-injection current. 

The critical current suppression by spin injection has 
been observed in FM/SC heterostructures made of a 
high-T c SC and a ferromagnetic manganite with P ~ 
100% @@. The experimental result of I c is surpris- 
ingly similar to that shown in Fig. 4. This strongly 
suggests that the injection currents from FM build up 
the spin density in SC of the heterostructures. Since 
the spin density is accumulated more efficiently in the 
double-junction geometry, the FM/SC/FM double junc- 
tions using high-T c SCs and ferromagnetic manganites is 
quite promising to test our predictions. 

Another candidate for SC in the FM/SC/FM junc- 
tion is a thin film of clean SC with sufficiently long 
spin-relaxation time such as Al. The depression of the 
TMR effect is caused predominantly by spin relaxation 
due to spin-orbit scattering in SC. Here, we derive a 
condition for observing the TMR effect when SC is in 
the normal state. By balancing the population rate 
(III — /2|)/ e with the relaxation rate AcW^Sa/ts, where 
d is the thickness of SC and A the junction area, we ob- 
tain Sa = ^PeS^V with the effective spin polarization 
Pes = P/(1 + r t /r.) where r t = AdV N h(R FN /R K ) || 
with R K = h/e 2 w 26 kft and i? FN = 1/G FN - To re- 
tain a substantial value of P e g, we are required to satisfy 
Tt < t s . For the case of clean Al with the values of 
t s - 10~ 8 s § and V N ~ 10 22 /(eV-cm 3 ), we have the 
condition Rp^A < (63/d) x 10~ 5 SI cm 2 for the specific 
junction resistance i?FN^4 and d (A). 

In the above calculations we have not considered the 
effect of Andreev reflection (AR) on the tunnel conduc- 
tance, since the resistance of a tunnel junction with a thin 
insulating layer is much higher than that of a metallic 
contact, i.e., Rp^A 3> 7ri?K:/2fc§- ks 10~ n ft cm 2 , where 
hp is the Fermi momentum. In metallic contacts with 
resistance comparable to 7rRjt/2kp, on the other hand, 
the conductance is dominated by AR for eV < 2 A [ pi] . 
Recently, the suppression of AR in FM/SC nanocontacts 



has been used to measure P in ferromagnets pj§ . 

In summary, we have studied the spin-dependent tun- 
neling in a FM/SC/FM double tunnel junction. The 
spin imbalance in the tunneling currents gives rise to the 
nonequilibrium spin density in SC. The superconductiv- 
ity is strongly suppressed with increase of bias voltage 
and destroyed at the critical voltage V c (oc 1/-P). The 
tunnel magnetoresistance exhibits an unusual voltage de- 
pendence around V c below T c . The results predicted in 
this Letter provide a method for measuring the spin po- 
larization of FMs as well as for controlling superconduc- 
tivity and TMR by application of bias voltage. 
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